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ABSTRACT 

Context. A population of obscured supergiant High Mass X-ray Binaries (sgHMXBs) has been discovered by INTEGRAL. X-ray 
wind tomography of IGR Jl 7252-36 16 inferred a slow wind velocity to account for the enhanced obscuration. 
Aims. The main goal of this study is to understand under which conditions high obscuration could occur 

Methods. We have used an hydrodynamical code to simulate the flow of the stellar wind around the neutron star. A grid of simula- 
tions was used to study the dependency of the absorbing column density and of the X-ray light-curves on the model parameters. A 
comparison between the simulation results and the observations of IGR J17252-3616 provides an estimate on these parameters. 
Results. We have constrained the wind terminal velocity to 500-600 km s~^ and the neutron star mass to 1.75 - 2.15 M©. 
Conclusions. We have confirmed that the initial hypothesis of a slow wind velocity with a moderate mass loss rate is valid. The mass 
of the neutron star can be constrained by studying its impact on the accretion flow. 

Key words. X-rays: binaries. Hydrodynamics, Stars: winds, outows. Accretion, accretion disks. Stars: individual: IGR /17252-3616 



1. Introduction 

In classical supergiant High Mass X-ray Binaries (sgHMXBs) 
neutron stars are orbiting at a distance of a ~ 1.5 - 2 R* from 
their companion stars. The donors, in these systems, are OB su- 
pergiants with mass loss rates of the order of ~ 10"^ M© yr"^ 
and wind terminal velocities of ~ 1500 km s"^ . 

[Blondin et al.| ( |1990l|1991| ) modeled the interactions between 
the neutron star and the stellar wind in Vela X-1 and revealed 
that the wind of the massive star is disrupted by the gravity and 
photoionization of the neutron star. The imprint of these two pa- 
rameters can be observed in the variability of the absorbing col- 
umn density with orbital phase. High-resolution soft X-ray spec- 
troscopy of the brightest sources revealed a number of lines in 
emission, constraining the ionization level of the gas ( jWatanabe 
let al.|2QQ6| ). 



The heavily obscured sgHMXBs ( Walter et a l. 2006) share 
some of the characteristics of the classical sgHMXBs. The main 
diflTerence between classical and obscured sgHMXBs is that the 
latter ones are much more absorbed in the X-rays. The absorbing 
column density (Nh ^ 10^^ cm"^) is, on average, 10 times larger 
than in classical systems and well above the galactic absorption 
in the direction of the sources. 

The obscure d sgHMX B IGR /172 52 - 3616 (=EXO 1722- 
3 63; [Walter et al. 2006; Zu rita Heras et a l. 2006 : Thompson et al. 
|2Q0 7') is an eclipsing binary hosting a pulsar with -414 sec, 
an orbital period of P^, ~ 9.74 days, and an orbital radius of 
a ^1.15 R^. Ground-based observations ( |Mason et al.|20Q9a|b| ) 
showed that the donor star is likely a B0-5I or BO-1 la. Optical 
and IR observations confir med the supergiant nature and showed 
prominent P-Cygni profile ( [Chaty et al.|2008| ) 



To explain the X MM - Newto n and INTEGRAL observa- 
tions, |Manousaki£_&_^^lter hereafter MWll) suggested 
that the wind terminal velocity of the system is relatively, low, of 
the order of Voo ~ 500 km s"^ An ad-hoc modeling of a hydro- 
dynamic trailing tail allowed to reproduce the observed column 
density profile and the observed iron Fe Ka line emissivity. 

In the present paper, we have used the hydrodynamic code 
VH-1 to study the wind dynamics of IGR J17252-3616 under 
the assumption of a slow wind velocity. We start with a descrip- 
tion of the hydrodynamic code and of the grid of simulations in 
section [2| the results are presented in section [3] and discussed in 
section |4j followed by the conclusions in section |5] 



2. Hydrodynamical Simulations 

2.1. The hydrodynamical code 

The full description of the VH-l[^ hydrodynamical code is given 
Blondin et ah] ( [T99Q| [T99T] ). To simulate sgHMXB the code 



m 

accounts for: i) the gravity of the primary and of the neutron star, 
ii) the radiative acceleration of the stellar wind of the primary 
star, and iii) the suppression of the stellar wind acceleration in 
the Stromgren sphere of the neutron star. The simulations take 
place in the orbital plane, reducing the problem to two dimen- 
si ons. The code uses the piecewise parabolic method developed 
by lColella & Woodward] 84). 

A computational grid of 600 radial by 247 angular zones, 
extending from 1 to ~ 15 R=, and in angle from -tt to +7r, has 
been employed. The grid structure is non-uniform, with cells 



http :// wonka. phy sic s . nc su . edu/pub/ VH- 1 / 
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size decreasing towards the neutron star. The maximal resolu- 2.3. Photoionization 
tion, reached close to the neutron star, is SR ~ 10^^ cm ~ 



where rgcc = 2GMx/(v^^^ + v^) is the accretion radius dBondi & 



Hoyle|1944| ). 



The neutron star in IGR J 17252-36 16 orbits in a almost cir- 



cular (e < 0.15) and highly inclined (/ > 80°) orbit (Manousakis 
& Walter] |2011| ). Therefore, a circular (e = 0) and edge-on 
(/ = 90°) orbit is assumed for simplicity. As most of the ab- 
sorption comes from the accretion stream, close to the neutron 
star (within 10^^ cm), the variability of the absorbing column 
density with orbital phase does not depend on small variations 
of the inclination angle. 

The code produces density and ionization = Lx/nr^^, 
where Lx is the average X-ray luminosity, n is the number den- 
sity at the distance from the neutron star; Tarter et al.|1969| ) 
maps that are stored. These allow to determine the simulated 
column density. As short time-scale variations occurs, we have 
calculated the time-averaged orbital phase resolved column den- 
sity and its 3cr variability. The 3-D instantaneous mass accretion 
rate (Mace) onto the neutron star is also recorded by translating 
the 2-D instantaneous mass accretion rate along the rim of the 
accretion radius, race- From the mass accretion rate, we can in- 
fer the instantaneous X-ray luminosity of the neutron star. For 
each simulations, the code was run for about 6 to 8 orbits, corre- 
sponding to about 60 - 80 days. This is enough for the wind to 
reach a stable configuration. The relaxation time is of the order 
of 0.8 orbits. The first 10 days of the simulations are therefore 
excluded from our analysis. 

Although 3D hydrodynamic simulations are poorly studied 
in HMXB, diff'erences between 2D and 3D are expected. In 
general, 3D simulations increase the instabilities and fluctua- 
tions. These efl'ects have been studied in the idealized case of 
Hoyle-Lyttleton accretion ( |Taam & FryxeU|1988[|Ruff'ert|1999[ 
Blondin & Raymer|2012| ) and in accreting pulsars where a for- 
mation of a planar accretion disk is expected from 2D codes, 
whilst accretion in 3D is unstable, allowing the formation of fil- 
amentary structures ( [Ruff'ert & Amett|1994 ). 



2.2. Stellar wind acceleration 

The winds of hot massive stars are characterized, observation- 
ally, by the wind terminal velocity and mass-loss rate. The ve- 
locity is described by the yS- velocity law, v = Voo(l - R^/r)^ 
where Voo is the terminal velocity and /3 is the gradient of the 
velocity field. For supergiant stars, values for wind terminal ve- 
locities and mass-loss rates are in the range Voo ~ 500 - 3000 km 
s~^ and Mw ~ 10"^ - 10"^ M© yr"\ respectively (Kudritzki & 
Puls|200 ). The donor star typically weights > 10 M© with a 
radius of ^ 10^^ cm. 

The winds of massive supergiant stars ar e radiatively driven 
by absorbing photons from the photosphere ( Castor et al.|1975 



CAK hereafter). The modeling of stellar wind is complicatec 
by the fact the rad iation force is known to be unstable ( Qwocki 
& RybicEl|1984 ). We use the CAK/Sobolev approximations 
with the finite disk correction ( [Friend & Abbott||1986| ) which 
produces a smooth stellar wind with JS = 0.8. However, re- 
gions in the stellar wind can diff'er from the predictions of the 
CAK/Sobolev approxima tion when gr owth of instabilities are 
taken into account ( |Owocki et al.|1988"] ). 



The radiative acceleration of the wind from the donor star pri- 
marily occurs due to the UV line transitions that are suppressed 
by X-ray ionization. The eff'ects of X-ray ionization on the radia- 
tive acceleration force are complicated due to the large number 
of ions and line transitions which contribute to the UV opac- 
ity ( [Stevens & Kallmaii||1990| ). A critical ionization parameter 
can be defined, above which the radiative force is negligible. For 
^ > 10^-^ erg cm sec"\ most of the elements responsible for the 
wind acceleration are fully ionized and, therefore, the radiatively 
acceleration force vanishes. At that level, light atoms (e.g. H) 
are fully ionized and do not contribute to the absorbing column 
density while heavier atoms (e.g. C, N, O) are highly ionized 
( [Kallman & McCray|1982| ) 

The main eff'ect of the ionization is the enhancement of the 
mass accretion rate onto the compact object. The acceleration 
cutoff' also triggers t he formation of a dense w ake at the rim of 
the Stromgren zone ( [Fransson & Fabian||1980 ), and has an im- 
pact on the absorption at late orbital phases. X-ray ionization can 
also aff'ect the thermal state of the wind through X-ray heating 
and radiative cooling, which can create additional filamentary 
structures in the wind, with slightly diff'erent time-averaged ab- 
sorbing column density ( [Blondin et al.|199Q l. 



2.4. Parametrization 

In order to simulate IGR J17252-3616, we proceed in two steps. 
We started by investigating the parametrization of the stellar 
wind of the donor star. This is achieved using a 1-D simu- 
lation assuming spherical symmetry. Spectral classification al- 
lowed to constrain the donor star parameters ( [Chaty et al.|2008[ 
Mason et al.[[2009a|b] ). In our analysis we have adopted an ef- 
fective temperature of Te/f = 30000 K, a stellar luminosity of 
L* = 4 X 10^ Lq, and a stellar radius of R^ = 29 Rq. The mass 
of the donor star is assumed to be M^, = 15 M© ( |Takeuchi et slI] 
T99Q[ [Thompson et al.pOOT] ). 



The wind acceleration depends on the CAK-^ and CAK-of 
parameters representing the number and strength of the absorp- 
tion lines, and po, the density at the surface of the donor star. 
Each of the triplet (k, a, po) produces a monotonically increas- 
ing wind velocity. Validation of the wind is performed by fitting 
the yS- velocity law to the simulated data which returns the wind 
terminal velocity and the gradient of the velocity field. The j3 
parameter was found to be ~ 0.8 in all cases and the terminal ve- 
locity ranged between 500 and 1200 km s"^ The corresponding 
mass-loss rate (using the continuity equation) is ranging from 
2 X 10"^ Mo yr-i to 5 x 10"^ M© yr'^ Table [l] shows the CAK 
parameters and the corresponding mass-loss rates and terminal 
velocities. Although the stellar parameters are fixed in our sim- 
ulations, observational uncertainties are present. The simulated 
absorbing column density profile can be aff'ected by these uncer- 
tainties (see sect. |4.2| for more details). 

We used the above CAK parameters to trigger the 2-D sim- 
ulations. The system now hosts a neutron star of mass Mx (in 
Mq), having a luminosity Lx = 10^^ erg s"^ (MWll), placed at 
a distance a (in R^) from its donor. A grid of parameters with 
varying a and was used to simulate the system. The param- 
eters are listed in table [T] 
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Table 1: Parameters of the 2-D simulations. Three different groups are considered: Top: Comparison between slow and normal wind 
terminal velocity for a number of mass-loss rates and different neutron star masses. Middle: Simulations to study the dependency, 
of the most likely model, on the neutron star mass. Bottom: Fine tuning of the neutron star mass and binary separation. 



Model name^ 


group 


CAK-Qf^ 


CAK-^^ 




b 

Po 






d 






vl2_ML2_A175_MNXYZ 




0.55 


0.05 


1 


10" 


-12 


2 


1200 


1.75 


1.55, 1.75, 1.95 


vl2_ML10_A175_MNXYZ 


1 


0.55 


0.15 


5 


10- 


-12 


10 


1200 


1.75 


1.55, 1.75, 1.95 


V 1 2_ML 5 0_A 1 7 5_MNXYZ 




0.55 


0.30 


1 


10" 


-11 


50 


1200 


1.75 


1.55, 1.75, 1.95 


v5 . 0_ML2_A175_MNXYZ 




0.35 


0.20 


5 


10" 


-12 


2 


500 


1.75 


1.55, 1.75, 1.95 


v5 . 0_ML10_A175_MNXYZ 




0.35 


0.30 


5 


10" 


-11 


10 


500 


1.75 


1.55, 1.75, 1.95 


v5 . 0_ML50_A175_MNXYZ 




0.35 


0.55 


1 


10" 


-10 


50 


500 


1.75 


1.55, 1.75, 1.95 


v5 . 0_ML10_A175_MNXYZ 


2 


0.35 


0.30 


5 


10" 


-11 


10 


500 


1.75 


1.5, 1.6, 1.7, 1.8, 1.9, 2.0 


v5 . 5_ML10_A173_MNXYZ 


3 


0.37 


0.28 


5 


10" 


-11 


10 


550 


1.73 


1.85, 1.90, 1.95, 2.00 


v5 . 5_ML10_A174_MNXYZ 


3 


0.37 


0.28 


5 


10" 


-11 


10 


550 


1.74 


1.85, 1.90, 1.95, 2.00 


v5.5_ML10_A175_]y[NXYZ 


3 


0.37 


0.28 


5 


10" 


-11 


10 


550 


1.75 


1.85, 1.90, 1.95,2.00 


v5 . 5_ML10_A176_MNXYZ 


3 


0.37 


0.28 


5 


10" 




10 


550 


1.76 


1.85, 1.90, 1.95, 2.00 


v5 . 5_ML10_A177_MNXYZ 


3 


0.37 


0.28 


5 


10" 




10 


550 


1.77 


1.85, 1.90, 1.95, 2.00 


V 6 . 0_ML 1 0_ A 173 .MNX YZ 


3 


0.38 


0.27 


3 


10" 


-11 


10 


600 


1.73 


1.85, 1.90, 1.95, 2.00 


v6 . 0_ML10_A174_MNXYZ 


3 


0.38 


0.27 


3 


10" 


-11 


10 


600 


1.74 


1.85, 1.90, 1.95,2.00 


v6 . 0_ML 10_A 1 7 5_MNXYZ 


3 


0.38 


0.27 


3 


10" 


-11 


10 


600 


1.75 


1.85, 1.90, 1.95,2.00 


v6 . 0_ML10_A176_MNXYZ 


3 


0.38 


0.27 


3 


10" 




10 


600 


1.76 


1.85, 1.90, 1.95, 2.00 


v6 . 0_ML10_A177_MNXYZ 


3 


0.38 


0.27 


3 


10" 




10 


600 


1.77 


1.85, 1.90, 1.95, 2.00 


v6. 5_ML10_A173_MNXYZ 


3 


0.395 


0.25 


3 


10" 


-11 


10 


650 


1.73 


1.85, 1.90, 1.95, 2.00 


v6. 5_ML10_A174_]y[NXYZ 


3 


0.395 


0.25 


3 


10" 


-11 


10 


650 


1.74 


1.85, 1.90, 1.95,2.00 


vfi S ML10 A17S MNXYZ 


3 


0.395 


0.25 


3 


10" 


-11 


10 


650 


1.75 


185 1 90 1 95 2 00 

X.O.^, X.VVy, X.V.^, Z^.WVy 


v6. 5_ML10_A176_MNXYZ 


3 


0.395 


0.25 


3 


10" 


-11 


10 


650 


1.76 


1.85, 1.90, 1.95,2.00 


v6. 5_ML10_A177_MNXYZ 


3 


0.395 


0.25 


3 


10" 


-11 


10 


650 


1.77 


1.85, 1.90, 1.95,2.00 


v7 . 0_ML10_A173_MNXYZ 


3 


0.417 


0.22 


2 


10" 


-11 


10 


700 


1.73 


1.85, 1.90, 1.95,2.00 


v7 . 0_ML10_A174_MNXYZ 


3 


0.417 


0.22 


2 


10" 


-11 


10 


700 


1.74 


1.85, 1.90, 1.95,2.00 


v7 . 0_ML10_A175_]y[NXYZ 


3 


0.417 


0.22 


2 


10" 


-11 


10 


700 


1.75 


1.85, 1.90, 1.95,2.00 


v7 . 0_ML10_A176_MNXYZ 


3 


0.417 


0.22 


2 


10" 


-11 


10 


700 


1.76 


1.85, 1.90, 1.95,2.00 


v7 . 0_ML10_A177_MNXYZ 


3 


0.417 


0.22 


2 


10" 


-11 


10 


700 


1.77 


1.85, 1.90, 1.95,2.00 



(a) Model name used as a reference; XYZ= 100^ Mx. (b) The triplet of input parameters CAK-of, CAK-k, and density, po, at the surface of the 
donor star in gr cm~^. (c) The mass-loss rate in 10~^ M© yr~^ (d) The wind terminal velocity in km s~^ (e) The binary separation in R^. (f) The 
masses of the neutron star in M©. 



3. Results 

3.1. Wind terminal velocity and mass-loss rate 

In order to study and constrain the dependency of the absorbing 
column density on the wind properties, we built a grid of simula- 
tion for a variety of mass-loss rates and wind terminal velocities. 
We have used two wind terminal velocities (500 km s"^ and 1200 
km s"^). and mass-loss rates of (0.2, 1, 5) xlO"^ M© yr"^. The 
orbital radius is fixed to a=l.l5 and the mass of the neutron 
star varies in the range 1.55, 1.75, and 1.95 M©. Six wind config- 
urations (CAK-Qf, CAK-k, Po) were simulated for each of these 
masses. These parameters are listed in table [T] (group 1) 

The geometry of the shock, formed around the neutron star, 
is affected by the wind parameters and by the ionization. The 
position and orientation of the bow shock depends on the relative 
velocity of the neutron star within the wind and on the X-ray 
luminosity. 



Figure [T] shows the density maps (after 3 orbits) for some 
wind configurations. The orientation of the bow shock and ac- 
cretion wake (tail) depends on the mass-loss rate and wind ter- 
minal velocity. For the slow wind, the accretion wake is always 
visible, enhanced, and tilted. For a wind velocity of 1200 km 
s"\ the accretion wake is not that tilted. For a mass loss rate of 
M = 2 X 10"^ Mq yr"\ the accretion wake is much weakened. 

Figure [2] shows the time-averaged simulated absorbing col- 
umn density as a function of orbital phase, for Mx=1.95 M© 
with the observations of IGR /17252- 3616. All these profiles 
show some common characteristics. At early phases, the smooth 
wind component dominates, followed by a rise of Nh (orbital 
phase ~ 0.3), representing the bow shock and the accretion wake. 
Then, the absorption roughly remains at a constant level. The 
tidal stream between the donor and neutron star also influences 
the absorbing column density at later phases. 
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Among these curves, only one simulation can represent the 
data (v5 . 0_ML10_A175_MN195) suggesting a wind terminal ve- 
locity of Voo = 500 km s"^ and a mass loss rate of M = 10"^ M© 
yr"^ The dependency of the absorbing column density orbital 
profile with the mass of the neutron star and binary separation 
parameters are discussed in sections [T2| and [33) respectively. 
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Fig. 3: Time-averaged simulated absorbing column density (Nh) 
for various neutron star masses (group 2 in table [T]). The stellar 
wind is characterised by a M = 1 x 10"^ M© yr"^and 1^00=500 
km s"^ The data points are from MWl 1. 



3.3. Binary separation 



Fig. 2: Time-averaged simulated absorbing column density (Nh) 
(group 1 in table [T]). The mass of the neutron star is fixed to 
Mx=1.95 Mq. The data points are from MWll. 



3.2. Mass of the neutron star 

A group of simulations was implemented with neutron star 
masses varying between 1.5 and 2.0 Mq. The summary of 
these simulations is given in table [T] (group 2). The separation 
and the wind parameters are identical to those of the model 
v5.0_ML10_A175_MN195. 

Figure [4] shows the density maps (after ~ 3 orbits). The cor- 
responding time-averaged absorbing column density is plotted 
in figure [3] The heavier the neutron star, the more bended the 
accretion wake. This eff'ect is revealed in the orbital phase de- 
pendency of the absorbing column density, moving the position 
of the minimum to earlier phases. For an heavier neutron star, 
more gas will accumulate and the absorption will get stronger. 
We can therefore estimate, in the frame of our model, that the 
mass of the neutron star in in the range M^s = 1.9 - 2.0 M© (see 

fig.©. 



The wind structure is also aff'ected by the binary separation. 
When the neutron star is closer to the donor star, the tidal stream 
between the two stars is enhanced. This eff'ect is maximized 
when the donor star fills its Roche lobe. In classical and obscured 
sgHMBs it is very likely that the evolved companion is close to 
filling its Roche lobe. Therefore, tidal connectio n between the 
two stars enhances the absorbing column density ( Blondin et al. 
fT99Tl). 



Figure [5] shows the absorbing column density for various 
wind terminal velocities . In these plots there is an indication 
that for a combination of neutron star mass and binary separa- 
tion we sustain a roughly constant obscuration. Table[2]gives the 

of the simulation (from fig. [5]) accounting for both the error 
bars of the data and temporal variations of the absorbing col- 
umn density. The data can be represented by the model with a 
wind terminal velocity, neutron star mass, and binary separation 
in the range of 500-600 km s-\ 1.85-1.95 M©, and 1.74-1.75 R,, 
respectively. 



3.4. Statistical analysis of the X-ray light-curve 

We have compared the observed X-ray light-curves of IGR 
/1 7252 - 3616 with these derived from our simulations. The 
observed light-cur ves of IGR J 17252 - 3616 were obtained form 
the INTEGRAL ( |Winkler et al.|2003| ) soft y-ray imager ISGRI 
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density (gr/cm-^) 




density Cgr/cmA3) 



X (xlO"!! cm) 




X (xWlZ cm) 



(a) Vo, 



500 km s-i ; ~ 2 X 10"^ Mq yr" 



(b) z^oo ^ 1200 km s-^ ; M^. ~ 2 x 10"^ M© yr" 



density (gr/cmA3) 



-3. 

I 




density Cgr/cmA3) 
-l.ODOe-13 




(c) Voo ^ 500 km s-i ; ~ 1 X 10"^ Mq yx'^ 



(d) Voo ^ 1200 km s"^ ; ~ 1 x 10"^ M© yr" 



Fig. 1: Density distribution (gr cm ^; color bar) in the plane of the orbit after ~ 3 orbits for a neutron star mass of Mx = 1.95 M© 
and different wind parameters. 



Table 2: The;^^ for models for various neutron star masses, bi- 
nary separations, and wind terminal velocities. 



Model name 



X 



v5 


O_ML10_A176_]y[N200 


3.5 


v5 


0_ML10_A175_MN195 


0.9 


v5 


0_ML10_A175_MN190 


1.8 


v5 


0_ML10_A175_MN185 


1.9 


v5 


0_ML10_A174_MN195 


1.3 


v5 


0_ML10_A174_MN190 


0.8 


v5 


0_ML10_A174_MN185 


1.1 


v5 


0_ML10_A173_MN195 


1.5 


v5 


0_ML10_A173_MN190 


2.3 


v5 


0_ML10_A173_MN185 


2.8 


v5 


5_ML10_A175_MN195 


0.8 


v6 


0_ML10_A175_MN195 


1.7 


v6 


5_ML10_A175_MN195 


3.0 


v5 


5_ML10_A174_MN190 


0.7 


v6 


0_ML10_A174_MN190 


1.3 


v6 


5_ML10_A174_MN190 


1.9 


v5 


5_ML10_A173_MN185 


1.5 


v6 


0_ML10_A173_MN185 


1.4 
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( [Lebrun et al.||2003] ) i n the energy range 20-60 keV using 
HEAVENS0interface ( [Walter et al.|2oTo| ). 



Light-curves were built using data from the source and from 
a source free near-by region, using 1-hour, 3-hour, 6-hour and 
12-hour time-bins and flux histograms were constructed. To 
achieve a relatively narrow width for the background fluctuations 
and a good time resolution, we selected the 6-hour time-bins for 
our analysis. 

The histogram of the observed fluxes can be fitted with nor- 
mal (gaussian) distributions. The observations can be character- 
ized by a mean count rate of yu=l.l cps and a width cr = 1.0 
cps for the source and yu=0 cps with a width cr = 0.6 cps for the 
background (black and green lines in fig.|6]). 

We performed a statistical analysis of the simulated light- 
curves which have been built using the mass accretion rates de- 
rived from the hydrodynamical simulations, translated into the 
instantaneous luminosity using Lace = Tl^accC^^ where // = 0.1. 
One should notice that the input parameter Lx - 10^^ erg 
s"^ = const, (see sect. [2]) and the instantaneous X-ray luminosity 
are not the same. However, the average Lace is within the 10% 
of the input parameter Lx- The Lace histogram can be fitted by 
a log-normal distribution (see inset of fig. [6]). The lognormal fit 
resulted in a luminosity ~ 1.1 x 10^^ erg s"^ having a width 
of cr = 0.44 for model v5.5_ML10_A175_MN195. 



\protect\vrule width0pt\protect\href {http : //www 



To compare the simulated distributions with the observed 
ones, we convolved the simulated data with the background fluc- 

^j^g tuations (blue curves in fig. [6]). We, also, plotted the simulated 
Nh for this model together with its 99% temporal variability 
(right panel of fig.|6|). 

isdc . unige . ch/heavensHhttp : //www . isdc . unige . ch/heavens} 



5 



A. Manousakis et al.: Neutron star masses from hydrodynamical effects in obscured sgHMXBs 



density Cgr/cm'^) 



I 




X (xlO'^12 cm) 



--I 4.0- 




X (xl0^12 cm) 



(a) Mx = 1.5 Mo 



(b)Mx = 1.7Mo 



density (gr/cmA3) 
].000e-12 



-2. 

t 




(c) Mx = 1.9Me 



X (xl0^12 cm) 



(d) Mx = 2.OM0 



Fig. 4: Density distribution (gr cm ^; color bar) in the plane of the orbit after ~ 3 orbits with a terminal velocity of Voo ~ 500 km s ^ 
and a mass-loss rate M^^ ~ 10"^ M© yr"^ The mass of the neutron star scales from 1.5 to 2.0 M©. See the text for more details. 
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Fig. 5: Time-averaged absorbing column density (Nh) for various parameters (see group 3 in table[T]). The stellar wind is character- 
ized by M = 1 X 10"^ M© yr"^ The data points are from MWl 1. 



4. Discussion 

4.1. Wind terminal velocities in obscured systems 

The wind terminal velocities of supergiant OB stars is corre- 
lated with the effective temperature (for a review see|Kudritzki & 
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Cmnt Hate (cpt^l 



Fig. 6: Results from the simulation |v 5 . 5_MLlQ_A175_]y[N195| Left: The histogram of observed light-curve (black), the background 
fluctuation (green), and the convolution of the simulated data with the latter (blue). The inset shows the corresponding distribution 
of the simulated, light-curve (black) with a log-normal fit (red dashed curve) . Right: The corresponding simulated time-averaged 
absorbing column density, together with its 99 % temporal variability, and the observed data from MWl 1. 



Pu1s|20QQ"] ). For supergiant stars wind terminal velocities scales 
from 300 km s"^ (for T^ff =10 kK) to 3000 km s'^ (for Teff =40 
kK). An imprint of these outflows is observed in the P-Cygni 
profiles of UV/Optical lines. For the system IGR J17252-3616, 
Qptical/IR studies a re diflScult due to high interstellar absorption 
dRahoui et al.l2008l ). 



For IGR J17252-3616, we suggested a wind terminal veloc- 
ity of 500-600 km s"^ which is about 2-3 times less than ex- 
pected for an OB supergiant of the same spectral type and ef- 
fective temperature. The wind in massive binaries can be slower 
than in isolated OB stars. The neutron star can indeed cut-off' 



the acceleratio n by ionizing the local environment ( [Stevens & 
Kallman|199Ql ). Note that the pecuHar X-ray binary, GX 301-2, 
hosting a massive early type star (Wray 977), also features a low 
wind terminal velocity of Voo ~ 500 km s"^ ( Kaper et al.|2006 ). 



The column density in this system reaches ~ 10^"^ cm"^ and is 
interpreted a s the signature of a high density stream ( [Leahy &] 
Kostka|2 ()08^). In this context, it is also interesting to remark that 
some intermediately obscured sgHMXBs, such as 4U 1907 -h 09, 
feature intermediate wind terni inal velocity in the range ~ 1000 
km s-i ( [Kostka & Leahy|2010t . 



4.2. The mass of the neutron star 

The neutron star mass determined by the hydrodynamical simu- 
lations is independent from the mass function derived from dy- 
namical studies. However, it is model-dependent and uncertain- 
ties on the modelization of the wind may have an impact. Figures 
[3] and [5] show that the mass of the neutron star is in the range 
Mx = 1.85 - 1.95 Mq. This mass is below the theoretical max- 



imum mass of a neutron star, M^^^ < 3 M© ([Miiller & Serot 



[T9961 ). 

The mass of the neutron star and the orbital radius are de- 
generate as seen in fig. [5] However, this degeneracy cannot allow 
very large or small masses of the neutron star (see table |2]). The 
eclipse duration constrains the binary separation from 1.7 to 1.8 
R:, (MWll) and therefore the mass can scale from 1.85 to 2.0 
Mq. a rather heavy neutron star is needed in any case to account 
for the observed profile. 

Radial velocity studies of the source, using VLT, provided 
estimates on the masses of the neutron star and donor star in the 
range Ma^5 = 1.46+0.38 M© and Moe = 13.6+1.6 M© for Roche 
lobe overffow; M^s = 1.63 + 0.38 M© and Mqb = 15.2 + 1.6 M© 



for an edge-on view ( [Mason et al.|2009b| ). In both cases the mass 
of the neutron star is consistent with our results 
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The companion radius and mass are obtained from the du- 
ration of the ecHpse and stellar classification from mid-IR spec- 
troscopy. These values are subject to observational uncertainties. 
To study the eff'ects of these uncertainties on the determination 
of the mass of neutron star, we ran simulations of our best model 
for M, = 15 ± 1 Mo and =29 + 3 Rq. The eff'ect of the stellar 
radius on M^s scales as M^siR*) ^ 0.0697?*. Given the range 
of acceptable stellar radius (i.e., 26 to 32), the resulting neutron 
star's mass range varies from 1.75 M© to 2.15 M©. The impact 
of the donor's mass on Ma^^ scales as M^si^*) ^ 0.13 M*, re- 
sulting in, Mns = 1.8-2.1 M©. 

Other supergiant HMXBs also feature large neutron star 
masses. Vela X-1 host s a neutron star with mass in the range 
Mx = 1.88 ± 0.13 Mo dQuaintrell et al.|2003|) and 4U 1700-302 
in the range Mx = 2.44 + 0.27 Mq ( |Clark et al.|2002| ). Although 
the sample is small, masses of neutron star in sgHMXBs, favors 
normal nucleonic EOSs ( |Lattimer & Prakash|2007| ). 
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INTEGRAL Workshop "The Restless Gamma-ray Universe" (INTEGRAL 

2010), 162. 

Walter, R., Zurita Heras, J., Bassani, L., et al. 2006, A&A, 453, 133 
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5. Conclusions 

In this paper we compare the observed properties with the results 
of hydrodynamic simulations of the heavily obscured sgHMXB 
IGR J17252-3616. Our conclusions are: 

- A low wind terminal velocity, Voo = 500 - 600 km s"\ is 
needed to explain the observations. A higher wind terminal 
velocity cannot produce enough obscuration to account for 
the observations. 

- The neutron star mass can be constrained from the variability 
of the absorbing column density along the orbital phase. In 
the case of IGR J 17252-36 16 the mass of the neutron star is 
in the range 1.75 - 2.15 Mq. 

- The separation of the system can be limited in the range 
a=1.73-1.75 R,. 

Measuring the variability of the absorbing column density 
with the orbital phase in obscured sgHMXBs provides an inde- 
pendent constraint on the mass of neutron stars. 
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